Relic Abundance Predicts Universal 
Mass- Width Relations for Dark Matter Interactions 
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We find new and universal relations for the properties of dark matter particles consistent with 
standard relic abundances. Analysis is based on first characterizing the s-channel resonant anni- 
hilation process in great detail, keeping track of all velocity-dependence, the presence of multiple 
scales and treating each physical regime above, below, and close to thresholds separately. The res- 
onant regime as well as extension to include non-resonant processes are then reduced to analytic 
formulas and inequalities that describe the full range of multi-dimensional numerical work. These 
results eliminate the need to recompute relic abundance model by model, and reduce calculations 
to verifying certain scale and parameter combinations are consistent. Remarkably simple formulas 
describe the relation between the total width of an s-channel intermediate particle, the masses and 
the couplings involved. Eliminating the width in terms of the mass produces new consistency rela- 
tions between dark matter masses and the intermediate masses. The formulas are general enough 
to test directly whether new particles can be identified as dark matter. Resonance mass and total 
width are quantities directly observable at accelerators such as the LHC, and will be sufficient to 
establish whether new discoveries are consistent with the cosmological bounds on dark matter. 

PACS numbers: 95.35.-|-d 



Thermal evolution of dark matter in the early universe 
can be found in textbooks to predict a velocity-averaged 
annihilation cross section (cru) 3 x 10~^^cm'^/s. How- 
ever the textbook exercise happens to assume velocity- 
independent cross sections that are neither general, nor 
reliable. In this paper we give a complete and general 
analysis based on a new strategy. Our method solves the 
inverse problem of bounding the multi-dimensional pa- 
rameter regions such that the relic abundance is fixed. 
The analysis involves several novel steps. 

Many models are summarized by t and u-channel ex- 
changes that are slowly varying, plus s-channel reso- 
nances that are a great complication. An early study 
by Greist and Seckel [l| noted that resonant processes 
violate the assumptions of constant cross sections, while 
being impossible to summarize with equivalently simple 
formulas. A typical resonant calculation involves several 
coupling constants and 5 dimensionful scales: the incom- 
ing energy, two masses mx, my, the final state mass or 
masses, and the final relic density. Choosing one point 
in this huge parameter space and solving the Boltzmann 
relic evolution will predict one particular relic density. 
Complete exploration has previously appeared impossi- 
ble, and most studies are limited to checking consistency 
in a plane of a few selected parameters. 

Our approach first characterizes relic evolution via s- 
channel annihilation. The work is partly numerical, and 
partly analytic; we make no assumptions about masses, 
couplings, or final states. We find several tricks to iden- 
tify the important scales and ratios of scales that describe 
every possible parameter region. We use unitarity and 
the optical theorem to represent exactly the s-channel 
decay into all possible final states. 



As a result, we find universal mass- width relations 
which fit the numerical work across the whole param- 
eter range. Let X be the dark matter (mass mx), and 
Y be the s-channel intermediate connector (mass my, 
width Ty), by which X + X ^ Y ^ anything. The 
mass-v-width relation for the pole below threshold is 
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Here axxY is the coupling of dark matter to the inter- 
mediate connector particle, and Cjj' is a spin-counting 
coefficient from which dimensionful scales have been re- 
moved. An equivalent formula for a pole above threshold 
{my > 2mx) is presented in Eq. [161 Section ITVl 

To show how the relation works, suppose an s-channel 
connector of mass my — 600 GeV has coupling axYY = 
10-2 and width Ty = 5 GeV. Then mx = 417.8 GeV 
is the mass of the dark matter if the pole lies below 
threshold and dominates the relic evolution. This is far 
more precise and specific than the order-of-magnitude es- 
timates generated by (av) 3 x 10~^^cm^ /s. If we have 
the same mass mx , and the width Ty = 2 GeV, then the 
resonance Y is too narrow to give the usual relic abun- 
dance. 

Many examples of dominant s-channel dark matter an- 
nihilation models can be found in the literature [2t |3|- 
However we are not limited to s-channel dominance. The 
most general cross section, including any number of chan- 
nels and interference, is either larger or smaller than the 
s-channcl annihilation. Supposing that the cross section 
with extra non-resonant channels is larger provides an 
inequality of a particular sense, given in Section [II B I An 
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inequality of the opposite sense comes from the oppo- 
site assumption. In general the s-channel dominant case 
produces a bounding surface in the space of all the pa- 
rameters. The surface is simple enough that a number of 
powerful inequalities in selected parameter planes come 
out rather easily. But we can do more: Eq. |9] shows 
quantitatively how to take into account any amount of 
resonant versus non-resonant cross sections with a simple 
"replacement rule." 

The mass-width relations also serve as a test of any 
new physics compared to dark matter cosmology. If the 
width and mass of a resonance do not match our relations 
then it is not a candidate to produce relics. Conversely a 
match of mass and width would be an indisputable signal 
of discovery. Note the width of a new particle is always a 
physical observable available from its production, allow- 
ing direct data-versus-data comparisons not depending 
on model details. Thus we expect our formulas to be 
useful for dark matter studies at the LHC and the ILC 



Moreover, widths are calculable in almost every model. 
The ansatz Ty = army is typical, but we also use it as a 
definition of the symbol ar for analysis. Recall that the 
Standard Model Z-boson has a mass of mz = 91.2GeV 
and a total width Tz = 2.1GeV, giving an effective cou- 
pling of az = 0.03. The effective coupling absorbs aU the 
channels, including the invisible ones, all the couplings, 
spin factors, and phase space. Combining Ty — array 
with Eq. [T] yields a non- linear equation relating nix-, rny 
and the couplings. For mx >> 100 GeV the relation 
reduces to my = 2mx + ci + 02/ mx + ■■■ where ci, C2 
are known functions of couplings. It is very surprising 
there is always an allowed solution for arbitrarily large 
mx, my, contrary to the expectations of Born-level esti- 
mates. In the event that symbol ar contains some mass 
dependence, as with certain theories with dimensionful 
couplings, the mass relation remains good and can be 
explored without needing to repeat the relic abundance 
calculations. 

We built on previous experience with s-channel anni- 
hilation effects in the galactic halo ff]. Velocity depen- 
dence in annihilation came to the front with recent satel- 
lite data from PAMELAJI^ERMI Q, PPB-BETS [3 
and other experiments ll,[ll|. Whether or not the data 



might be a signal, the studies have led to recognition that 
Born-level cross sections are not adequate. Exaggerated 
claims about "Sommerfeld factors," sometimes thought 
to be exact non-perturbative effects, both violate general 
principles and phenomenological tests [l3l - IT5j . Refs. 
01 flfil [l7| found the enhancements of ordinary resonant 
physics can be surprisingly large. Resonant processes can 
saturate unitarity bounds on annihilation in halo circum- 
stances, generating large "boost factors" suggested by the 
data. It is because certain width-dependent effects ,per- 
turbatively small in the high energy limit, may dominate 
everything in the non-relativistic limit. The problems 
of hevy thremal relic evolution in the early universe are 
similar, because the abundance is determined primarily 



in a regime sensitive to non-relativistic dynamics. Refs. 
(isl [l9| show concrete examples in which even the effects 
of one loop corrections on the relic abundance can be 
significant. 

There are many applications of our results. Section |T] 
reviews the standard relic abundance formalism for com- 
pleteness. Section |ll] presents the s-channel mass- width 
relation for a pole below threshold, along with motiva- 
tion for the analytic formula. At each stage where it 
is appropriate we convert the equalities developed for 
s-channel dominance into inequalities, illustrated with 
graphics. Section Hill explores consequences of calculable 
widths. This leads to mx — my relations between the in- 
termediate particle and the dark matter. The surprising 
fact that my 2mx can fit relic abundance for arbitrar- 
ily large masses is something to explore. It turns out to 
be a generalization of the "funnel' region" known in mini- 
mal supersymmetric models (SUSY) [20l - [23 . [26j . without 
making an assumption of SUSY. We suspect that our re- 
lations are more restrictive than the commonly known 
ones, and may rule out some models. Separate analysis 
for the case where the intermediate connector mass my 
is above the kinematic threshold 2mx is given in Section 

llYl 



CROSS SECTIONS AND RELIC 
ABUNDANCE 






FIG. 1: s-channel dark matter annihilation diagrams into all 
possible final states. 



The thermally evolving number density Y is calculated 
as a function of the inverse temperature x = mx /T using 

dY 

— mx){av){x){Y'-YiQ), 

where Yeq is the equilibrium density, and ^{x, mx) = 
—xs{x)/H{mx) is a combination of standard entropy 
and Hubble functions [131 ■ The asymptotic solution fol- 
lows in the regime of Yeq 0: 

dY 



— = ^(a-, mx){(Jv){x)dx; 



dx^{x, mx){o'v). 



(1) 



The lower limit Xd is computed self-consistently. Given 
the standard value of the critical density, the asymptotic 
density VL = mxYooSo/ pcrit follows, using the critical 
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density Pcrit and entropy density sq of the present uni- 
verse. 

Particle physics enters in the annihilation cross section 
a{v) and its thermal average {(jv){x). The rapid energy 
dependence and multiple scales of resonant cross sections 
require special analysis. Expansions of the form {av) k, 
a + 6/x + c/x^ + ... are never good approximations over 
the entire range of x. 

Let Mxx^f be the amplitude for XX to go to a final 
state /. The cross section a goes like the amplitude- 
squared, summed over all final states (FigH]), and inte- 
grated over final state phase space LIPS: 
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The total cross section into all possible final states is 
given by the optical theorem: 
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Here k is the momentum of either particle in the center of 
mass frame, and M is the elastic scattering amplitude. 
For a given total center of mass energy Ecm and its 
square s, the forward propagators of intermediate states 
Y go like (s — m|, + iTOyr)~^, where F is the total width. 
Let gxxY^jj' be the component of the elastic amplitude 
containing the couplings of the initial/final states of spin 
j to an s-channel particle of spin j' . Then channel by 
channel, the optical theorem predicts 
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The symbol Fy represents the total width of Y to all final 
states, which allows us to describe numerous models with 
a single parameter. 

In standard convention for amplitudes, the Feynman 
rules contain in/out state polarization and vertex factors 
compiled into the symbol tjji . It is important to extract 
the mass (mx) dependence of these factors for analysis. 
We define 



X'-'jj'^ 



(3) 



whereby Cjji is typically a number of order unity. Scaling 
like is expected from dimensional analysis, and in- 
evitable when the initial state is dominated by the mass 
as the largest scale. We emphasize that Cjj' is a defi- 
nition that allows for any model while postponing spin- 
sums and vertex factors until a model is chosen. For 
example, the annihilation of unpolarized Dirac Fermions 
via a 7'' vertex produces Ci/21 — 3/4. Another exam- 
ple is the annihilation of two scalar particles to a vector, 
which will have two derivatives going like k^ ^ m\ , times 
the polarization sum over the vector particle modes. The 
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FIG. 2: Intersection of the mesh plane representing Q.h^ =0.1 
and the surface r2/i^(ry, my, m,, a*) gives a unique curve 
r(my). m, = IQQGeV , q, = 0.01 for the purpose of the 
graphic. 



coupling and spin factors then appear in the combination 
A-KaxxyCjji , where axYY = q^xyH'^'^)- 

Now all mx dependence of the cross section, rate {crv), 
and density Yoo dependence has been scaled out, except 
for a minor cutoff {xd) dependence in Eq. [1] which must 
be reconciled numerically. In the numerical work we 
first choose a particular dark matter mass mx = ti* 
and coupling a*, and then compute thousands of relic 
densities covering the Fy, my plane. The condition 
ri(Fy, m*, my, a*)/i^ — > 0.1 produces a unique curve Fy 
versus my, namely the function Fy — Fy(my; m*,a*), 
as seen in Fig. [2] With the scaling relations in hand, 
the curves are extended to numerical predictions for the 
general functional dependence of F = F(my; mx, otxyy) 
consistent with a fixed relic density. 



II. MASS- WIDTH RELATIONS: POLE BELOW 
THRESHOLD (my < 2mx) 

Fig. [3] shows the mass- width relation as a family of 
curves plotted for selected mx- The trend is that the 
further the my is from the threshold, the larger Fy must 
be needed to keep relic densities constant, and vice versa. 
This is because the proximity to threshold (rather than 
the absolute size) of my is the dominant effect. Poles 
closer to the threshold make for larger cross sections, 
which need to be compensated by smaller width. 

This quantitative understanding leads to remark- 
able analytic formulas reproducing the whole parameter 
range. 



A. Analytic Representation 

Observe in Fig. [3] that each dashed (black) curve 
moving to the right terminates in a region of Fy — )■ 0. 
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mY(GeV) 

FIG. 3: Relation of the s-channel width Fy and mass my 
for a pole below threshold consistent with cosmological rehc 
density Qh^ =0.1. Dashed curves (black online) are the nu- 
merical calculation. Solid curves (blue online) are the analytic 
relation of Eq. |6] Each curve is evaluated with a fixed dark 
matter mass mx =100-500 GeV in 100 GeV increments. 



Two terms suffice with gi = 2. Our analytic formula for 
the pole below threshold mass-width relation is then 
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The fit of the formula to the numerical work is extremely 
good. Fig. [3] shows a typical example. 

We have already noted that Eq. IH] generalizes and re- 
places the traditional formula (av) = 3 x 10~^^cm^ /s. 
The formula accounts for the fact that intermediate 
states of all particles coupling to dark matter are either 
absolutely stable or have a finite lifetime. To underscore 
the difference, compare the traditional counting rules, 
motivated on dimensional analysis, using annihilation 
cross sections of order 1-picobarn. Under the assumption 
av ^ /'^x ^ P^i ^ typical upper limit a\ < lO"** 
would imply mx ^ 200 GeV. This well-known result 
needs a finely tuned mass-couping relation to make a 
Universe. Our formula contains far more information, 
and reveals the hidden assumption that my « mx was 
implicitly assumed for the traditional formula to be con- 
sistent. 



Near the threshold everything is determined by the de- 
gree of the zero of the function, Ty ~ [my — 2mx)". 
The power n — 2 can be gotten analytically, but was 
also fit directly with numerical work. Then we know 
Fy (1 — my /2mxY , times a known factor of mx 

In the opposite extreme of my << 2mx the veloc- 
ity averaged cross section {av) reduces to another simple 
analytic result: 
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Notice the dependence going like l/{myTy) 
the equation gives 
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The limit of small connector mass my — >■ with 
Fy << 2mx approaches the Born approximation, and 
for us is the unique case where the Born cross sec- 
tion is relevant. In the Born limit we know {av) 
3 X 10~^^cTO'^/s, which fixes one overall scale. Then ac- 
counting for factors of mx gives 
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The dimensionless interpolating function g{my /mx) re- 
mains. It must obey g — 1 when my << 2mx, suggest- 
ing a polynomial expansion g ^ 1 + X^fe 9k{'^Y I'n^x)^ ■ 



Replacement Rule for Adding Non-Resonant 
Channels 



Some models of dark matter annihilation include more 
than one channel. No matter how many channels are 
involved, as long as the s-channel is a part of the model, 
the mass-width relation of Eq. [5] can be generalized. 

For definiteness, suppose the addition of other pro- 
cesses increases the annihilation rate. Then the theory 
keeping Q.h'^ fixed will require smaller Fy, all other things 
fixed. That condition rules out all contours to the right 
and above the contours shown in Fig. [3] created by s- 
channel dominance. The allowed region to the right and 
below each line implies an inequality (Fig. ^ : 
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While cross sections often increase when channels are 
added, destructive interference occurs in some models. 
In that case the inequality reverses the sign. 

To illustrate the use of Eq. [71 suppose a new vector 
boson {Z' perhaps) is discovered at the LHC. Measur- 
ing the resonance observed in any channel will give its 
mass my, and the total width Fy. Applying Eq. [7] then 
gives the consistent mx, ctxxY parameter space regions 
consistent with relic abundance without a need for an ex- 
tensive numerical parameter space scan. These predicted 
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The revised mass-v-width relation then becomes 



200 300 400 500 600 700 800 
mx(Gev) 

FIG. 4: Upper limits on Vy assuming s-channel annihilation 
(pole below threshold) plus other channels increasing the cross 
section. Shaded regions to the left and above the contours 
are not allowed. Curves show different couplings axxY = 
5xxv/47r; my = 400 GeV is used for the purpose of the 
graphic. Larger my pushes contours to the right. 



coupling relations are then compared to the information 
from production rate and branching ratio into particular 
channels seen in the experiment. 

A new relation comes from a "replacement rule." Let 
the total velocity averaged cross section be expressed as 



the 



Suppose (o'w)Qj^gj. happens to be consistent with the tra- 
ditional Born-style of approximation, by which 
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Fig. [5] shows the replacement rule performs quite well. 
Not surprisingly, the difference relative to a pure s- 
channel annihilation model increases for smaller mx- 
This is because the individual cross section contributions 
from other channels shown in Eq. [9] scale uniformly like 
\lra\. 




200 400 600 800 1000 
nriY (GeV ) 

FIG. 5: Generalization of the mass-width relation to include 
"Born-like" channels. Black dashed curves show numerical 
evaluation. Solid curves (magenta online) are the revised fit of 
Eq. [9] Thick solid curves (blue online) are the approximations 
of Eq. [51 consistent with the role as an upper bound. Different 
curves use different masses mx = 100 — 500GeV^, from left to 
right. Parameter ctxxY = 10~^ and X^i^te// — 10~^ for the 
purpose of the graphic. 



where a^^j is an effective coupling to the i^^ channel. 
(An example model in the context of heavy hidden sec- 
tor dark matter can be found in Ref. [23].) Matching 
the extreme limits produces the replacement rule. If the 
s-channel pole is near the threshold, the mass width re- 
lation should approach Eq. [B] If the pole is far from 
threshold the resonant cross section approaches an effec- 
tive Born- level cross section, and adds to it. That implies 
a boundary condition of {av) ^ -I- {<yv)^^f^^^ « lO'^Gel^"^ 
at this endpoint. Reviewing how that scale previously 
entered the analysis suggests a replacement rule: 

2.6 X 10-9Gey-2 ^ 2.6 x lO-^Gel/-^ _ (g) 



III. CALCULABLE WIDTHS CONSTRAIN THE 
MASSES 

Up to here we have considered the width Ty as an 
independent parameter. In this section we go a step fur- 
ther and consider widths as quantities which can be cal- 
culated. When we say that "widths are calculable" it 
emphasizes the facts that (1) most theories are pertur- 
batively coupled, and (2) most of the width will usu- 
ally occur in a finite number of channels. Whatever the 
models combining the calculation of the width with the 
mass- width relation creates a new relation. 
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The differential rate dT of a general decay of a particle 
of mass TOy is given by 



dT 



2mY 



IMl'^dLIPS. 



Symbol M is the amplitude. The final state phase space 
of two identical particles yield / dLIPS2 ^ Vf where Vf 
is the velocity of either final state particle in the center 
of mass frame. In many cases the width is dominated by 
relativistic final states, w/ — J- 1. It would be unusual, and 
a case of rather fine tuning, for all channels with phase 
space limitations Vf << 1 to dominate the total width. 
Barring that event, by dimensional analysis, the width 
of a heavy particle with dimensionless coupling generally 
goes like its mass: 



Fy ~ ariTiY- 



(10) 



We make this an equality allowing symbol ar to ab- 
sorb coupling constants, the number of important chan- 
nels, and model details. The general scaling of widths- 
proportional-to-mass is rather kinematic. However, if 
a dimensionful coupling is introduced, then the mass 
dependence of rest of the calculation is dominated by 
dimensional analysis again. Keeping in mind that ar 
stands for the width actually calculated in a particular 
model, we continue. 

With ar fixed, the formula for Fy is an increasing 
function of my. Meanwhile the s-channel mass- v- width 
requirements are all decreasing functions of my (Fig. |3l 
Eq. ini Eq. [5]). Then the width- v- mass relation always 
matches the s-channel mass-v-width relation at a definite 
point. Fig. [6] shows Fy = apwy as red curves, whose 
intersections with the blue curves constrain the masses. 

Inspection finds a surprising fact. Rather weakly cou- 
pled theories (ar < 10^'^) only intersect the relic curves 
in the region where my w 2m x- In minimal supersym- 
metry, this result corresponds to the so called higgs " fun- 
nel" region of mo,mi/2 parameter space. However, our 
result is much more general and extends beyond the as- 
sumptions of SUSY. 

The two couplings, ar and axxY act in the same di- 
rection. Smaller ar makes smaller widths that force the 
system into the threshold region to be viable. Smaller 
axYY worsens the situation by pushing the contours of 
constant f2/i^ up in the Fy — my plane. The trend of both 
pushes masses into very near coincidence of my ^ 2mx , 
which we call a "finely-tuned threshold." 

Except for bound state formation we have no reason 
to consider finely-tuned thresholds very plausible, but 
we can afford to stay neutral. Bound states have been 
discussed in detail in Refs. [l^lsO]- Some basic relations 
between bound state widths and masses are reviewed in 
Ref. 01 • Bound state relations are very specific and 
require separate treatment that is not our topic here. 

Our mass-width relation allows for classification of 
models according to the degree of fine tuning. For ex- 
ample. Fig. [5] shows a theory with axxY = 0.1. The 
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200 400 600 800 1000 
my (GeV) 

FIG. 6: Combining the below-threshold width-v-mass my 
relation of Fig. [3] with Ty = army represented by solid 
thin curves (red online). Intersections of the curves predict 
a non-linear relation between mx and my (text). Values of 
10"^ (top curve) range to ar = 10 * (bottom curve) in 
factor of 10 increments. 



intersections are not very demanding, and the theory is 
not finely tuned for ar > 10"'^. However the same theory 
using axxY = will require widths 100 times bigger 
for the same my to keep O/i^ constant. At that point 
all contours are pushed up to such a degree we'd find the 
theory finely-tuned for all reasonable Fy. 



A. Dark Matter and Pole Mass Relations 

So far we have looked at an my - Fy relationship, 
given mx- It is very interesting to consider the rela- 
tion between my and mx given ar. Fig. [7] shows the 
mass relationships for diff'erent values of ar. Once again 
ar forces the resonance into the finely tuned region 
of TOy « 2mx- 

Solving Eq. [6] yields a cubic equation fixing my = 
myimx , axxY , ctr)- The relationship is nearly linear 
for a wide range of mx ■ Collect the couplings into a new 
symbol 

2 OixxYCtrCjj' 



10-4 

Note the symbol has been re-scaled in units of 
ctxxY /W^^ , ar/lQ-^ we find reasonable. The series ex- 
pansion for large mx is found to be 

my « 2m,x - 3l3GeVa^ + m.2GeV al ( ) 

mx 

.5.87GeValC-^^r + ...in) 

mx 
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FIG. 7: Mass of dark matter mx versus the mass of the 
particle in the s-channel my- Red lines represent ar ~ 
0.001,0.01,0.1 from top to bottom. axxY = 0.01 for the 
purpose of the graphic. Small widths (ar small) require fine 
mass tuning, my ~ 2mx to accommodate correct relic abun- 
dance. 



Eq. [TT] with only the first two terms kept is essentially 
exact for a„ < 1, mx <; 100 GeV, while for > 1 a 
numerical evaluation is preferable. 

Analyzing Fig. [7] and [8] we notice that the mass range 
of 100-500 GeV does not require extreme fine tuning for a 
reasonable range of perturbative couplings. On the other 
hand the regime of mx » lOOGeT^ seems to require a 
pole tuned very finely according to Eq. 1111 
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FIG. 8: Same as Fig. [7] but for an extended range of mx- 
Red lines represent ar ~ 0.001, 0.01, 0.1 from top to bottom. 
axxY ~ 0.01 for the purpose of the graphic. 



IV. MASS-WIDTH RELATIONS: POLE ABOVE 
THRESHOLD {my > 2mx) 
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FIG. 9: Relation of the s-channel width Fy and pole mass 
my above threshold consistent with cosmological relic density 
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Dashed curves (black online) are the numerical 
calculation. Solid curves (blue online) are the analytic rela- 
tion of Eq. 1161 Each curve is evaluated with a fixed dark 
matter mass mx =100-500 GeV in 100 GeV increments. 



The relic abundance calculation for a pole above 
threshold is complicated by a saddle point in the inte- 
gration of {(Jv). To begin we again consider the extreme 
limits. For my >> 2mx the velocity averaged cross sec- 
tion reduces to 



{av){mY » 2mx) 
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xxvC'jj' 



mx 



Ty 



(12) 



By construction this limit reproduces the Born-level es- 
timate, {(Tv) « IG^^Gey^^. Introduce a dimcnsionless 
function h to describe other limits, expressed by 
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AnaxxvCjj' 
xh{my, mx, A, axxyCjj-). 

We have normalized /i — > 1 for my » 2mx by ab- 
sorbing the overall normalization into 6.4 x 10~9GeV~^ , 
suggesting the ansatz 



h = l- 



ymx 
my 



(13) 



where rj is the measure of the "offset" of my from the 
threshold 2mx- Unlike the case of pole below threshold, 
the saddle point causes to be a function of mx and a 
relic scale parameter, which we call A. 
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FIG. 10: Mass of dark matter mx vs. the mass of the s- 
channel particle my for a pole above threshold, my > 2mx 
Red lines represent av = 0, 0.01, 0.05, 0.1 from top to bottom. 
ctxxY = 0.01 for the pmpose of the graphic. Small widths 
(ar small) require very fine tuning of masses, my ~ r\mx to 
accommodate correct relic abundance. 



The extreme of F gives more information about 
the function h. The Breit-Wigner factor can be approx- 
imated as 



(s-m^)2 + (myr)' 



— >■ 7r(5(s 



Ly). 



(14) 



The Boltzmann equation is then solved analytically in 
terms of error functions, predicting h and r\ in this limit: 



r\{mx,oixxY,Cjji) 



—erfc i( ^ 




(15) 

where A — 1.3 x 10 ^^GeV ^ gives a good fit for all 
reasonable my and Fy. The lovi^er integration limit Xd 
is computed in a self consistent way and we find that the 
standard value of Xd = 30 is appropriate. 

Eq. 1151 involves the inverse complementary error func- 
tions (er/c~^), which is somewhat cumbersome. While 
many numerical packages (including Mathematica) com- 
pute it, a simpler analytic formulation of r] is useful. Let 
z = A /cexxyCjj' . We find the approximation 

r]{z) w 1.978 - 0.521Z - 0.051Log[z] 

is almost exact in the range 10~* < z < 1. 

Our analytic formula for a pole above threshold is now: 



V 



IGeV / my 
iTTaxxvCjj, \589GeVJ 

for my > 2mx 



1 - 



■n(z)mx 
my 



Once again, the analytic approximation matches numer- 
ical work remarkably well. Fig. [3] shows an example. 

Eq. [in] reveals more finely tuned parameter regions 
for a pole above threshold. In the limit Fy — >■ 0, my 
is finely tuned to rjuix, as seen in Fig. 1101 The com- 
petition between the pole position, width, and thermal 
Gaussian are all summarized by this generalization of the 
pole below-threshold relation. 

As before, eliminating Fy = army produces an mx ~ 
my relation: 



my 



r]{z)mx , 1 



2 2 
where ^ Sgg^GeF^ 



IQt: iJ? ar a xxyC 



'33' ' 




(16) 



200 300 400 500 600 700 800 
mx(GeV) 

FIG. 11: Typical upper limits on Fy given a mass my for 
different masses of dark matter mx assuming both s-channel 
annihilation (pole above threshold) and other non-resonant 
channels. Shaded regions to the left and above the contours 
are not allowed, my = 1.5 TeV was used for the purpose of 
the graphic; larger my pushes contours to the right. 



1. Upper Limit on mx 

The particular form of the "offset function" 77 yields an 
upper limit on mx- From the derivation the argument 
of erfc^^ must be less than one, which implies 

mx < 2.77 X lO^GeV^/axxvCjf (17) 

This formula is more precise than supposing my is 
bounded by a Born-level estimate and 2mx < my- 
Consider for example a small coupling axxy = 10^"*. 
Consistency with relic abundance requires dark matter 
masses mx ^ 2.8TeV. 
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2.. Inequalities for Non-Resonant Channels 

Generalization of the above-threshold mass width re- 
lation to allow for non-resonant channels is similar to the 
below-threshold case. When {av)^^^ > {<Jv) ^ the mass 
width relation becomes 



IGeV 



(- 



my 



ATraxxyCjf \589GeV 



(18) 



An illustration of the inequality can be seen in Fig. [TT] 
Notice that for large couplings, i.e. axxY > 0.1 major 
portions of the parameter space can be ruled out. The 
termination point {Ty ^ 0) is simply the rjinx = my 
point, giving us another bound on nix- By inspection, a 
coupling axxY =0.1 and my = 1.5TeV requires mx ^ 
600GeF. 



V. CONCLUSIONS 

The dynamical effects of resonant processes and finite 
particle widths play an important role in dark matter 



evolution in the early universe. The Born approximation 
is seldom adequate because the non-relativistic velocity 
dependence of cross sections drives decoupling. Orga- 
nizing the calculation in terms of observable quantities 
gives new relations between the masses and widths of in- 
termediate states that will be consistent with fixed relic 
abundance. 

Given that particle widths are generally calculable, our 
mass-v-width relations develop into mass-v-mass consis- 
tency relations between the dark matter with a given relic 
density and the mass of an s-channel connector. Depend- 
ing on the model, this produces a significant revision of 
a traditional rule (av) ~ 3 x 10~^^ cm^/s. The relation 
between mx and my depends on the way the width is 
calculated, but in a broad class of models permits an 
unlimited range of both masses. Our relations can be 
used to test candidates for dark matter in LHC-based 
experiments, while also eliminating much of the need to 
re-compute relic evolution on a model-by-model basis. 
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